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The integrals in Gradshteyn and Ryzhik. 
Part 11: The incomplete beta function 

Khristo N. Boyadzhiev, Luis A. Medina, and Victor H. Moll 

Abstract. The table of Gradshteyn and Rhyzik contains some integrals that 
can be expressed in terms of the incomplete beta function. We describe some 
elementary properties of this function and use them to check some formulas in 
the mentioned table. 



1. Introduction 

The table of integrals [?] contains a large variety of definite integrals that involve 
the incomplete beta function defined here by the integral 

f 1 r a ~ 1 dr 

(1.1) (3(a) = / X — ™. 

Jo l + x 

The convergence of the integral requires a > 0. Nielsen, who used this function 
extensively, attributed it to Stirling [?], page 17. The table [?] prefers to introduce 
first the digamma function 

(1.2) ^ ) = A logr(a; ) = EM ) 

and define P(x) by the identity 

(1.3) ^) = I( V) ( £ ±1)_ V ,(|)). 

This definition appears as 8.370 and 1)1. 1|) appears as 3.222.1. Here 

(1.4) r(a:) = / f-^dt 



is the classical gamma function. Naturally, both starting points for are equivalent, 
and Corollary 12.21 proves (|1.3p . The value 

(1.5) 7 := -V(l) = -I"(l) 
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is the well-known Euler's constant. 

In this paper we will prove elementary properties of this function and use them 
to evaluate some definite integrals in [?] . 

2. Some elementary properties 

The incomplete beta function admits a representation by series. 
Proposition 2.1. Let a e E + . Then 

■ (-If 



(2.1) P(a) = £ 



k=Q 



a + k 



Proof. The result follows from the expansion of 1/(1 + x) in (|1.1[) as a geometric 
series. □ 

Corollary 2.2. The incomplete beta function is given by 



(2-2) /3(a) = ~ 



This is 8.370 in [?]. 

Proof. The expansion for the digamma function ip 

°° f 1 1 
(2-3) lK*) = -7-£(F 

fc=0 

has been discussed in [?]. Then 



t+k k+1 



OO 

(2-4) ^(f)=-7-E 

k=0 

and 

(2-5) ^f^)=-7-£ 



2fc 



a + 2fc + 1 



The identity l|2.2p comes from adding these two expressions. □ 

These properties are now employed to prove some functional relations of the in- 
complete beta function. The proofs will employ the identities 

(2.6) ip(x + l) = -+ip{x) 

x 

(2.7) ip(x) - ip(l - x) = -7TC0t(7Tx) 

(2.8) tp(x + |) - ij)(\ - x) = Trtan(Trx) 

that were established in f?l . 
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Remark 2.1. Several of the evaluations presented here will employ the special 
values 

n 

(2.9) ifr(n + l) = -l + J2k> 

k=l 

that appears as 8.365.4, and 

(2.10) V (| ± n) = - 7 + 2 ^ ^ - ln2^ , 
that appears as 8.366.3. 

Many of the formulas in Section 4.271 employ the values 

2 n ~ 1 i 

that appear as 8.366.11 and also 8.366.12/13: 

2 n 1 

Higher order derivatives are given by 



V> (n) (l) = (-l)" +1 n!C(™ + l) and 
^(«)(|) = (-l)"+ 1 n!(2"+ 1 -l)C(n + l). 

Proposition 2.3. The incomplete beta function satisfies 

1 



(2.13) 
(2.14) 

(2.15) 

Proof. Using 
(3(x + 1) = 



f3{x + l) 
(3(1 -x) 

0(x + l) 

we have 
, fx + 2 



- P(x), 



sin 7tx 

1 7T 



/3(x), 



x sin 7rx 



+ j8(l-ar). 



-</> 



a: + 1 
2 

a; + 1 



03 + 1 



/3(x) 



This establishes ([2713"]) , To prove (|zTT4]l we start with 

/3(a:)+)9(l-i) - ~ 
The formula (|2.14|) now follows from (|2.7p and 



. 1 x 
^2 + 2 



□ 
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3. Some elementary changes of variables 

The class of integrals evaluated here are obtained from fll.l|) by some elementary 
manipulations. 

Example 3.1. The change x = t p in yields 

f 1 jtop-l 

(3.1) /3(a) =p ' 



Replace a by - to obtain 3.241.1: 



(3.2) /'^ = ;"(i 

Example 3.2. The special case p = 2 in Example 13.11 gives 

(3.3) « ) = 2 ^L_5. 

Choose a = 2±1 5 and relabel the variable of integration as x, to obtain 3.249.4: 

1 x b ax _l„fb + l 
o 1+x 2 ~ 2 

Example 3.3. The evaluation of 3.251.7: 



(3-4) / f^ = ^ 



[ x a dx 1 a— \ n (a— 1 

(3-5) „ : -^ = -T + — r-/? 



/o (1 + x 2 ) 2 4 4 r V 2 
comes from the change of variables t = x 2 and integration by parts. Indeed, 



x 



a dx 1 f 1 ,„ , w , d 1 



(1 + x 2 ) 2 2 J dtl + t 

1 O-l /"I f<*-3)/2 



t(»-i)/2 4_JL_dt 



4 4 ./n If/ 

and (|3.5p has been established. 

Example 3.4. Formula 3.231.2 states that 

( 3.6) r^ i+x - p dx =-^. 

J 1 + x sin Tip 

The integrals is recognized as (3(p)+(3(l—p) and its value follows from (|2.14[) . Similarly, 
3.231.4 is 

3.7 / — dx= ; . 

J 1 + x p sm irp 

The integral is now recognized as [3{1 + p) — /3(1 — p), and the result follows from 

(us}. 
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Example 3.5. The evaluation of 3.244.1: 

(j-p-i 



(3.8) 



is 



(3.9) 



1 + xi 



1 



7T piT 

dx = — cosec — 



~(J3{p/q)+p(l-p/q)) 



according to (|3.2p . The result now follows from (|2.14p . 
Example 3.6. The evaluation of 3.269.2: 



(3.10) 



x p — x 



-v 



■ dx 



1 



/ 1 + x 2 p 2sin(7rp/2) 

is obtained by the change of variables t = x 2 , that produces 



(3.11) 



1 



f p/2 _ t -p/2 



2 Jo 1 + t 
The result now follows from (|2. 151) . 



dt 



1 r 

2 







P 



1 



4. Some exponential integrals 

In this section we present some exponential integrals that may be evaluated in 
terms of the /3-function. 

Example 4.1. The change of variables x — e~* in (jl.ip gives 
(4.1) /3(a) = 



,o 1 + e 
This appears as 3.311.2 in [?]. 

Example 4.2. The evaluation of 3.311.13: 



(4.2) 



■ dx 



it / np 

-cosec 



/o l + e-(P+i) x p + q \P + q 

is achieved by the change of variables t = (p + q)x that produces 



1 



00 p -pt/(p+q) 



p + q Jo 1 + e ~ 



■dt- 



i 



00 p-qt/ip+i) 



1 



/5 



p + q 



+ j9 1 



P + q Jo 1 + 
P + g 



The result now comes from (|2.15p . 
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5. Some trigonometrical integrals 

In this section we present the evaluation of some trigonometric integrals using the 
/3-function. 

Example 5.1. The change of variables x = tan 2 t in gives 

pir/4 

(5.1) (3(a) = 2 1 tan 2 "" 1 tdt. 

Jo 

Introduce the new parameter b = 2a — 1 to obtain 3.622.2: 

W 4 1 fb -4-1 

(5.2) / t&n b tdt = -B 

Jo 2 

Example 5.2. The change of variables x — t&nt in (|3.5p gives 

/ \ r /4 „ ^ , i o-i /o-i 

(5.3) / tan Q t cos 2 tdt = 1 8 

Jo 4 4 

Now use (|2. 13|> to obtain 

(5.4) flf^ + ' 



a- 1 



that converts ()5.3p to 



(5.5) / tan a t cos 2 tdt = - -\ 8 

Jo 4 4 V 2 

This is the form in which 3.623.3 appears in [?]. Using this form and (|5.2p we obtain 
3.623.2: 

(5.6) r /4 tan a t s\n 2 tdt = -I + I±^/3 A ' + ' 
7o 4 4^ 

Example 5.3. The evaluation of 3.624.1: 

sin p xdx _ 1 

cos P+2 x p+l 

can be done by writing the integral as 



(5.7) 



ptt/4 / l7r /4 

(5.8) 1=1 t&n p+2 xdx + \ iwPxdx. 

Jo Jo 

These are evaluated using (|5.2j) to obtain 
The rule (|2.13p completes the proof. 
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Example 5.4. The integral 3.651.2 

(5.10) if, (!!+») t^"* 1 

Jo 1 — sin a; cos x 3 \ \ 2 / 

can be established directly using the integral definition of (3 given in Simply 
observe that dividing the numerator and denominator of the integrand by cos 2 x yields, 
after the change of variables t = tana;, the identity 

T/4 tan^xdx r /4 tan^x dx 



1— sin a; cos a; J (sec 2 x — tana;) cos 2 x 
A t" dt 



dt. 



t 2 -t + l 
i p+i + tf i 

t 3 + 1 

The change of variables t = s 1 / 3 gives the result. 
The evaluation of 3.651.1 

(5 id r /4 tsi ^ xdx = i (4 (t±i\ _ 

[ ' ' J Q 1 + sinxcosx 3 V \ 2 / V \ 2 

can be established along the same lines. This part employs the representation 8.361.7: 

(5.12) V(*)= f — = 1 - 7 

Jo as — i 

established in [?]. 

Example 5.5. The elementary identity 

(5.13) i = I ( I 

1 - sin a; cos 2 a; 2 \1 + sin a; cos a; 1 — smxcosx, 

and the evaluations given in Examples 15.111 and 15.101 gives a proof of 3.656.1: 

(5.14) £ - (^) + i, + ^ + 2* TO - 2^ (*±I)) 
Example 5.6. The final integral in this section is 3.635.1: 

(5.15) / cos Al_1 (2x) tanxdx = -/3(/x)- 
This is easy: start with 

sin x 2 sin x cos x sin 2x 

(5.16) tanx 



cos x 2 cos 2 x 1 + cos 2x ' 
and use the change of variables t = cos2x to produce the result. 
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6. Some hyperbolic integrals 



This section contains the evaluation of some hyperbolic integrals using the 
function. 



(6.1) (3(a) 



Example 6.1. The integral (|4.1|) can be written as 

/ e */2 + e -t/2' 
and with t = 2y and b = 2a — 1, we obtain 3.541.6: 

(6.2) f£^ =/3 ^±i 
V ' Jo coshy M V 2 

Example 6.2. Integration by parts produces 



2 / e- Qa: -?-- 



cosh 2 a; 7 dxl + e 2 ' 

= -l + 2a 



o 1 + e- 2x ' 

The change of variables t = 2x now gives the evaluation of 3.541.8: 
, x f°°e- ax dx /a\ 

Example 6.3. The change of variables t — e~ x gives 
(6.4) / e" aa; tanha;da; = / — — dt, 



and with s = t we get 



1 f 1 W 2 " 1 - s^- 1 )/ 2 

I = - ds 

2 in l+s 



1 

2 

The transformation rule (|2.13p gives the evaluation of 3.541.7: 

(6.5) / e- ax ta,nhxdx = (3(^) --. 

Jo \2J a 

7. Differentiation formulas 
Example 7.1. Differentiating (|l.ip with respect to the parameter a yields 

(7.1) f 1 ^ lnX dx = P'{a), 



1 + x 



that appears as 4.251.3 in [?]. 
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Example 7.2. Differentiating (|3.2[) n times with respect to the parameter a 
produces 4.271.16 written in the form 

Jo f + XP p n+L \p 

The choice n = 1 now gives formula 4.254.4 in [?]: 

f 1 x 11 " 1 \nx , 1 „, (a 

(7.3) / — dx = -zfi - 

J 1 + xP p 2 ' \ P/ 

Example 7.3. The special case n = 1 , a = 1 and p = 1 in (|7.2[) produces the 
elementary integral 4.231.1: 

1 In x dx 7r 2 



(7.6) / f^cfe 



(7-4) 

io 1+s 12 
In this evaluation we have employed the values 

(7.5) V'(l)=C(2) = y, and ^(1/2) = ^, 

that appear in (|2.12[) . 

Example 7.4. Formula 4.231.14: 

1 xlnx 
1 + x 2 "~ 48 

comes from (I7.3|) by choosing the parameters n = 1, a = 2 and p = 2. The values 
of ^'(1) an d i/j'(1/2) are employed again. Naturally, this evaluation also comes from 
(|7.4p via the change of variables a; 2 i— * x. 

Example 7.5. The choice n = a = 1 and p = 2 in (|7.3j) and the values 

(7.7) V (2) Q) - ^ 2 + 8G and V (2) (|) = vr 2 - 8G, 

where G is Catalan constant defined by 

yields the evaluation of 4.231.12: 

r In a; da; 
1 + x 2 

The change of variables x — t/a, with a > 0, and the elementary integral 

dt n 



(7-9) / -r—,=-G. 



(7.10) 

give the evaluation of 4.231.11 

(7.11) / 

Jo 



t 2 + a 2 4a' 
a Inxdx 7r In a — 4G 



x 2 + a 2 4a 
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Example 7.6. Now choose n = 1, a = 2 and p = 1 in (|7.3p and use the value 
V>'(3/2) = tt 2 /2 - 4 given in ([2~T2"1) to obtain 4.231.19: 

f 1 i lni , 7r 2 
(7.12) / dx = 1. 



Jo l+x 12 
Combining this with l|7.4p gives 4.231.20: 

f 1 1-X , , 7T 2 

(7.13) / lnxdx = l-— . 

17 Jo l + i 6 

Example 7.7. The values 

(7.14) V (2) (I) = -2tt 3 - 56C(3) and V (2) (§) = 2tt 3 - 56C(3), 

given in [?], are now used to produce the evaluation of 4.261.6: 

1 In 2 x dx 7r 3 
l + x ~ 16' 



(7.15) 

Example 7.8. The relation 



(7.16) ^ (n) (i - 2) + (-i)"+V (n) (^) = (-l)"*-— cot ^ z ' 



and the choice n = 4, a = 1 and p = 2 in ()7.3j) produces 

^ W (I) 



1 In 4 x dx 1 



1 + a; 2 2 5 

|i(> (4) (i)-^ (4) (l)) 



1024 
1 ( d 4 

— 7T— — r COt 7TZ 



1024 V dz 4 

This yields the evaluation of 4.263.2: 

j" 1 ln 4 xeb 57T 5 

(7.17) 

The evaluation of 4.265: 
(7.18) 



=3/4 



l + x 2 64 ' 

1 In 6 x dx 6l7r 7 



/ 1 + a; 2 256 ' 
can be checked by the same method. 

Example 7.9. Now choose n € N and take a = n + 1 and p = 1 in (|7.3|) to obtain 
the expression 

(7,9) ^'f^-^'^,,. 

This is now expressed in terms of the ^—function and then simplified employing the 
relation 

(7.20) ^ m \z) = (-l) m+1 m!C(m + l,z), 
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with the Hurwitz zeta function 

(7-21) «-,*):= Et^j. 

We conclude that 



The elementary identity 

(-l) fe 



(7-23) CM)-C(*,^)=2 s ]>: (fc _ r 

is now used with s = 3 and a = n + 1 to obtain 

f 1 x n \n 2 xdx _ n ^ (-l) k 
( • ' J l + x ^(fc + n+l)3- 

This is finally transformed to the form 

This is 4.261.11 of [?]. 

The same method produces 4.262.4: 

(7 26) ^^.^^..g^. 



Example 7.10. The method of the previous example yields the value of 4.262.1: 

1 In 3 x dx 77T 4 



(? - 27) ]„ 120- 

Here we use V (3) (!) = ^V 15 and ^ (3) (!/2) = tt 4 - 

Similarly, V (5) (l) = 8tt 6 /63 and ^ 5 \l/2) = 8tt 6 yields 4.264.1: 

f 1 In 5 a; da; 31tt 6 
(? - 28) 7 TT^ = - -252"' 

and V (7) (!) = 8tt 8 /15 and V (7) (l/2) = 136tt 8 yields 4.266.1: 

^Wxdx 1277T 8 

(7 - 29 » /„ — = ~535- 



12 
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(7.30) 

Indeed 



Example 7.11. A combination of the evaluations given above produces 4.261.2: 

1 In 2 x dx IOtt 3 



o 1-x + x 2 81a/3' 



1 ln 2 xdx f 1 1 + x 



In 2 x dx 



o 1 — x + x 2 J 1+x 3 

1 In 2 x dx f 1 x In 2 x dx 



o 1 + x 3 J a 1 + x 3 
^(/3 (2) (|) + /3 (2) (|)) 



d 2 

■ — — COt 7TZ 

216 \dz 2 

.2 



d 2 

+ -—^ COt 7T2I 

2=1/3 az^ 



2=1/6 



— — = + 8V3tt 2 



IOtt 3 



216 \Zy/Z / 81^3' 

Example 7.12. Replace n by 2n in ()7.2|) and set a = p = 1 to produce 
' ln 2n a;<2a; = ^(an)/^ 



2 2 

2 2n _ 2 



1 + x 

^(^(l)-^(l)) 
- 22 ^(2n)!C(2n + l). 
This appears as 4.271.1. 

Example 7.13. The change of variables f = bx in (|7.2|) produces 
6 /° -1 1n/ h a ~P / \ f b t a ~ l dt 



^'(f)+ &1 ~ a] 



bP + tP p 2 \PJ J bP + tP 

b a ~ b / \ 6 a ~ p / \ 



The last integral was evaluated using 

Differentiate this identity with respect to the parameter b to obtain 

(7 31) f *— lnt dt b0 ~ 2Plnb | p - a b«-*°3>(" 

J (bP + tP) 2M ~ 2p + p 3 b P \p 



P 2 



(1 + (a -p) hi &)/?(§) 
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The special case a = b = p = 1 yields 4.231.6: 
(T.32) 

Similarly, the choice a = 2, b = 1 and p = 2 yields 4.234.2: 

(7.33) / 1 S^ = -Z^) = - i ¥- 

Example 7.14. In this last example of this section we present an evaluation of 
4.234.1: 

In x dx G 7r 



(7 ' 34) < 2 

using the methods developed here. We begin with the change of variables x i— > 1/x to 
transform the problem to the interval [0, 1]. We have 

lnxdx f 1 x 2 \nxdx 



(7 ' 35) A (i + 7 (i + 

Now choose a = 3, 6 = — 1 and p = 2 in (|7.31j) to obtain 
fnn „s f 1 x 2 \nxdx l nl f3\ l„/3 

The value of (|7.34[) now follows from 

\HI) - IGKS) -*(§)) 

H4 + ^(|)-^(!)) 



and 



8 

1 _ 7T 

2 8 



8 



= ^(^(i)-^(!)-i6) 

= ^ (C(2, |) - C(2, |) - 16) 

_ G 1 

2 2 ' 

8. One last example 

In this section we discuss the evaluation of 3.522.4: 

f 00 dx _ 1 / 1 

' ' Jo (& 2 + x 2 ) cosher ~ b' \ + 2 



The technique illustrated here will be employed in a future publication to discuss many 
other evaluations. 

To establish (|8.ip . introduce the function 



oc 
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This function is harmonic and bounded for Re 6 > 0. Therefore it admits a Poisson 
representation 

If 00 b 

(8.3) h (b,y) = - H0,u)—— -du. 

7T J_ OQ b 2 + {y- u) 2 

The value h(0, u) is a well-known Fourier transform 

( 8 - 4 H0,u)= — ^-dt=- — — , 

y 1 Jo cosht 2 cosh(7r?//2)' 

that appears as 3.981.3 in [?]. Therefore we have 

hi* 00 rl 

(8.5) h(b,y) ' 1 



2 cosh(7ra/2) [b 2 + (y- u) 2 ] ' 
The special value y = and (|6.2|) give the result (after replacing b by 26 and u by 2u). 

Note 8.1. Formula (|8.ip can also be obtained by a direct contour integration. 
Details will be provided in a future publication. 

We conclude with an interpretation of (|8.ip in terms of the sine Fourier transform 
of a function related to (3(x). The proof is a simple application of the elementary 
identity 

f°° b 
(8.6) / e xt ambtdt = 



jo b 2 +x 2 ' 
The details are left to the reader. 

Theorem 8.2. Let 

, N [°° e~ xt dt fx+1 

^ ^--=1 ^hT^br- 

Then 3.522.4 in (|8.1|) is equivalent to the identity 



f°° /2x\ 
/ fi(t) sin xt dt = [i — I . 

V t J 
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